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We discuss the behavior of the extended sound modes of a dense binary hard-sphere mix- 
ture. In a dense simple hard-sphere fluid the Enskog theory predicts a gap in the sound 
propagation at large wave vectors. In a binary mixture the gap is only present for low 
concentrations of one of the two species. At intermediate concentrations sound modes 
are always propagating. This behavior is not affected by the mass difference of the two 
species, but it only depends on the packing fractions. The gap is absent when the packing 
fractions are comparable and the mixture structurally resembles a metallic glass. 



I. INTRODUCTION 

Short-wavelength collective modes in dense simple hard-sphere fluids have been stud- 
ied extensively by both theoretical and experimental researchers in the past few 
years. 1 " 4 Two interesting features of the extended hydrodynamic modes in a simple 
fluid are the softening of the heat mode at the wave vector k where the static struc- 
ture factor S(k) has its first maximum and the appearance of a gap in the sound 
propagation in the same large-wave-vector region. 1,2 The softening of the heat mode 
corresponds to the slowing down of structural relaxation in a dense fluid and has 
been discussed extensively in the literature. The sound-propagation gap consists 
in the vanishing of the imaginary part of the frequency of the sound waves, which 
describes the propagation of the wave, over a finite region of k values. 

In a one-component fluid of hard spheres of size a and density n one obtains a 
gap in the sound propagation for all values of the packing fraction 77 = (w/6)na 3 
above n m 0.33. 5 This feature of the extended modes has been the topic of some 
controversy, de Schepper, Cohen, and Zuilhof, 5 have shown that the appearance and 
disappearance of the sound propagation gap can be understood as arising from the 
competition between generalized wave-vector-dependent reversible elastic restoring 
forces and dissipative forces in the fluid. Roughly speaking, a gap arises whenever 
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the dissipation exceeds the elastic forces. They also pointed out that the behav- 
ior of the sound modes in a one- component fluid is quite insensitive to the static 
structure factor of the fluid and the gap is present even when setting S(k) = 1. 
de Schepper, Van Rijs, and Cohen 6 have also shown that even the long- wavelength 
Navier-Stokes equations, when used outside their range of validity to describe short- 
wavelength phenomena, predict a gap in the sound propagation. The model based 
on the Navier-Stokes equations suggests that there is no connection between local 
effects leading to wave-vector-dependent transport coefficients and susceptibilities 
and the trapping of sound on molecular length scales. 

In a recent paper (hereafter denoted I) we used the Enskog equation as the starting 
point to study the short-wavelength collective excitations of a binary hard-sphere 
fluid mixture. 7 In a mixture the spatial ordering of the particles can be changed 
while keeping the total density of the fluid constant by changing the concentra- 
tion X2 = 112/ {ni + n 2 ) of the larger species (here type-2 spheres) or the size ratio 
a = (J\ja-i- We find that sound propagation in the mixture depends on the con- 
centration x 2 : the trapping of sound modes and the corresponding propagation gap 
that arises in dense simple fluids only occurs in mixtures with a small concentration 
of small spheres in a dense fluid of large spheres (x 2 ~ 1) or a small concentration 
of large spheres in a dense fluid of small spheres (x 2 < 1). At intermediate values 
of X2 we observe a softening of the sound propagation, but no actual gap. For either 
small or large values of X2 the mixture closely resembles a simple fluid, as can be 
seen from the behavior of the partial structure factors Su(k), S22 (k) and S^ik) 
defined in I. For instance, for x 2 = 0.01 only the static structure factor Sn(k) of 
the dense component is peaked [see Fig. 1(a)]. The other structure factors are es- 
sentially constant. Thus we expect the appearance of a sound-propagation gap as 
in a simple fluid. At intermediate values of the concentration x 2 the dense mixture 
resembles a metallic glass: there is short-range order in the spatial arrangements of 
both types of spheres, as indicated by the fact that all three static structure factors 
display considerable structure as functions of the wave vector [see Fig. 1(b)]. The 
absence of a sound-propagation gap in this case is then consistent with the exper- 
imental finding that no gap occurs in the second sound of metallic glasses. 8 This 
observation seems to indicate that the local fluid structure does play a role in de- 
termining sound propagation and sound trapping at short wavelengths, in contrast 
to what was argued for a one-component fluid. 

In this paper we investigate in more detail the extended sound modes in binary 
mixtures, with the objective of clarifying their dependence on concentration, size 
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ratio and mass ratio, and the role of the fluid structure in determining the gap. 
As in our earlier paper, we keep the total packing fraction of the mixture rj = 
(7r/6)(ni0"f + n20"2) fixed at the value 77 = 0.46, which is close to the packing fraction 
corresponding to freezing of a one-component hard-sphere fluid. 

FIG.l. The partial static structure factors Sn(k), S22 (k), and ^(fc) for rj = 0.46, 
a = 0.7, x 2 = 0.01 (a), and x 2 = 0.5(b). 

First we have considered the dependence of the sound modes on the ratio mxjm^ of 
the masses of the two species and the role of this ratio in determining the propaga- 
tion gap. By comparing the extended hydrodynamic modes displayed in I [a = 0.7 
and mi/7712 = 0.5) to the modes obtained for the same values of size ratio and 
concentration, but equal masses, we have found that the sound propagation is es- 
sentially unchanged at all concentrations. This holds true for other values of the 
size ratio. In general, as long as the mass ratio is not too small (mi/7712 > 0.1) the 
difference in mass of the two species affects only weakly the sound propagation in 
a mixture. This points to the fact that the size difference and the spatial ordering 
of the two types of spheres must play the main role in determining the absence of 
sound trapping at large wave vectors. 

In Sec. II we discuss the dependence of the sound modes on the concentration of 
one of the two species. We show that in dense mixtures, as in a one-component 
fluid, the appearance of a sound propagation gap can be understood as the result 
of the competition between elastic and dissipative forces in the fluid. The absence 
of the sound-propagation gap in mixtures of intermediate concentrations is not de- 
termined by the details of the static structure factors, but it does depend on the 
fluid structure in the sense that it occurs when the packing fractions 771 = (7r/6)niof 
and t]2 = (7r/6)n2cr| of the two species are comparable. In this case all three static 
structure factors are peaked and the structural properties of the mixture resemble 
those of a metallic glass. 

II. DEPENDENCE OF SOUND PROPAGATION 
ON CONCENTRATION 

We showed in I that at large wave vectors the extended heat and diffusion modes 
govern the relaxation of the densities of the two species. These are the only long- 
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lived fluctuations in the fluid at these wave vectors. At large wave vectors the 
extended sound modes mainly describe the relaxation of temperature and longitudi- 
nal momentum fluctuations. Similarly, in a one-component fluid large-wave-vector 
density fluctuations are long-lived while temperature and momentum fluctuations 
relax quickly. On the basis of this observation Zuilhof and co-workerss suggested 
that extended sound modes in a one-component fluid can be described by a simple 
model obtained by neglecting the coupling to the density in the generalized hydro- 
dynamic equations for temperature and longitudinal momentum. The resulting two 
coupled equations have eigenfrequencies that closely reproduce the extended sound 
modes at large wave vectors. Clearly, at small wave vectors these modes do not 
reduce to the hydrodynamic sound modes. 

The same model can be used to describe the extended sound modes of a mixture, 
where the characteristic time scale for the relaxation of temperature and momentum 
fluctuations is well separated from that governing the relaxation of fluctuations in 
the two densities. If we neglect the coupling to fluctuations in the densities in the 
generalized hydrodynamic equations for temperature and longitudinal momentum, 
we obtain two coupled equations that are formally identical to those obtained for a 
one-component fluid. The eigenvalues of these equations are obtained by solving a 
quadratic equation, with the result 



where fizz (A;) and V^Trik) are the damping rates in the equations for longitudinal 
momentum and temperature fluctuations, respectively, and f(k) = iQTi(k) is the 
elastic restoring force that couples the two equations, 
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where n = n\ + n 2 and p = m\ni + 7712^2 are the total number and mass densi- 
ties, respectively. Also, p a b = , m a m b / {m a + nib) is the reduced mass, \ab the pair 
correlation function of species a and b at contact, a ab = {p- a + o"fc)/2, and j n (x) 
is a spherical Bessel function of order n. Finally, tE ab is the Enskog mean-free 
time, t Eab = y^JtabP / ' 4y/nn^nb(jl b Xab) ■ At large wave vectors the two modes given 
in Eq. (2.1) closely resemble the extended sound modes obtained in I by solv- 
ing the four coupled hydrodynamic equations. The argument of the square root 
in Eq. (2.1) is generally negative (and the two modes are propagating), unless 
the dissipative damping exceeds the elastic forces. In this case Eq. (2.1) yields 
two real roots, corresponding to diffusive modes. This is the region of the sound- 
propagation gap. The argument of the square root in Eq. (2.1) can be factorized as 
(flu — Qtt + 2/)(fi^ — Qtt — 2/). The first factor is always positive. We define the 
function A(k) by 



Q~12 

2 



3Pp 
2n 



1/2 



(n u -n TT -2f) = A(k)-k<T 12 , 



with 
A(fc) 



n 



n a n b a ab a 12 Xab 



0=1,2 6=1,2 



3np 



1/2 



j2(ka ab ) - jl(k<Tab) 



(2.5) 



A gap occurs when A(k) > koyi. This condition is displayed graphically in Fig. 2 at 
three different concentrations. The figure shows that a gap is obtained for x 2 = 0.01 
and x 2 = 0.9, but not for x 2 = 0.5. For a = 0.7 and mi/m 2 = 0.5 this simple 
model predicts no sound-propagation gap for 0.05 < 0.7, in agreement with what 
was obtained in I solving the full coupled hydrodynamic equations. 

FIG. 2. The straight line is ka\ 2 . The curves represent the function A(k) for 
77 = 0.46, a = 0.7, and m 1 /m 2 = 0.5, and x 2 = 0.01 (solid line), x 2 = 0.5 (dashed 
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line), and x<i = 0.9 (chain dashed line), as a function of ka 12 . 

In mixtures, as in a one-component fluid, the sound-propagation gap occurs as a 
consequence of the competition between dissipative and reversible forces. In most 
of the wave- vector range elastic forces exceed dissipative forces. This simply reflects 
the fact that sound modes propagate in both liquids and low-density gases (in the 
very-large- wave- vector limit the fluid resembles an ideal gas, since one is considering 
length scales that are too short for the interactions to be relevant). The damping 
term and the elastic forces depend on the packing fractions r\\ and 772. For a fixed 
total packing fraction rj and size ratio (a = 0.7), both dissipative and elastic terms 
decrease with increasing x 2 . At all x 2 there is a region of large k where the dissipa- 
tion rates are essentially constant (and equal the sum of the Enskog times weighted 
by the mass fractions) and Qu — Qtt ~ 0. For x 2 1 and £2 — the elastic forces 
vanish even more rapidly than the dissipative term in this large-fc region and a gap 
occurs. For intermediate x 2 the elastic forces remain finite at all k and the modes 
are always propagating. 

In a recent paper 9 Campa and Cohen suggested that the static structure factors 
may play a role in determining sound-propagation gaps in a mixture. On the other 
hand, the simple model described above does not contain the static structure fac- 
tors Sij(k). In fact, the appearance and disappearance of the sound-propagating gap 
are insensitive to the structure factors. This can also be seen by reconsidering the 
extended sound modes obtained in I from the solution of the four coupled hydrody- 
namic equation and arbitrarily setting Su(k) = S 22 (k) — 1 and Si 2 (k) = in these 
equations. One finds that, while the gap itself, when it occurs, is wider if the proper 
static structure factors are included, its concentration dependence, i.e., its appear- 
ance and disappearance, is qualitatively unchanged. The suggestion of Campa and 
Cohen is, however, in the right direction and it would be misleading to say that 
short-wavelength sound modes in a mixture are not affected by the fluid structure. 
The latter does enter even in the simple model yielding Eq. (2.1) through the de- 
pendence on the packing fractions of the two components of the mixture. When the 
packing fractions are both appreciable, there is not gap. This situation corresponds 
to a fluid where all three partial static structure factors are peaked (Fig. I) and the 
fluid structure resembles that of a metallic glass. 



III. CONCLUSION 
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We conclude with two remarks. 

(i) In I we neglected the coefficient of thermal diffusion in evaluating the extended 
hydrodynamic modes of the mixture. This approximation was motivated by 
the fact that the coefficient of thermal diffusion vanishes in a first-Sonine- 
polynomial approximation, and in the long-wavelength limit it is always much 
smaller than the diffusion coefficient in mixtures of spheres of not too disparate 
sizes and masses. 10 On the other hand, in the long- wavelength limit the coef- 
ficient of thermal diffusion contributes to the damping of the sound modes. 11 
Since neglecting thermal diffusion is best justified in mixtures of either very 
low 22 or 22 — 1, one may ask if the absence of the sound-propagation gap 
at intermediate 22 simply occurs because we have underestimated the sound 
damping by neglecting thermal diffusion. The model used here for the de- 
scription of the extended sound modes based on the two coupled equations for 
temperature and longitudinal momentum fluctuations neglects all couplings 
to the densities and therefore does not contain the coefficient of thermal dif- 
fusion. This model still predicts the disappearance and reappearance of the 
sound propagation g function of concentration, indicating that the de- 
pendence of the sound propagation on concentration is not an artifact of the 
approximation of neglecting thermal diffusion. 

(ii) Previous work has focussed on unusual sound propagation in disparate-mass 
gas mixtures and binary liquid alloys at moderate density. 9 ' 12 One interesting 
feature observed in computer simulations 13 is the appearance of a fast propa- 
gating sound mode above a certain nonzero value of the wave vector, signaling 
an effective separation of the dynamics of light and heavy particles. Campa 
and Cohen 9 have shown that the fast sound is a kinetic mode rather than a 
hydrodynamic mode and it yields an observable shoulder in the dynamic struc- 
ture factor of gas mixtures. In moderately dense mixtures the fast sound mode 
is present only for mass ratios below 0.1 and no fast sound occurs for rj > 0.42. 
Our generalized hydrodynamic equations as derived in I cannot yield a fast 
sound mode since our set of independent variables only includes the conserved 
densities of the fluid: we only obtain the extended hydrodynamic modes and 
no kinetic modes. To account for the possibility of fast sound we would need 
to enlarge our set of independent variables to include at least the relative mo- 
mentum density and the temperature difference of the two species. On the 
other hand, the fact that no fast sound is observed in Ref. 9 for the values of 
the total packing fraction and mass ratio considered here suggests that at such 
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high densities the dynamics of light and heavy particles never decouple and 
the generalized hydrodynamic equations of I may indeed provide an adequate 
description of the dynamics at large wave vectors. 14 A proper test of whether 
our theory is relevant for real mixtures will, however, only come from a de- 
tailed comparison of the dynamic structure factor that can be obtained from 
the generalized hydrodynamic equations of I with neutron-scattering spectra 
from dense mixtures or with numerical simulations. Finally, the fact that the 
appearance and disappearance of the sound-propagation gap in our work does 
not depend on the mass ratio and takes place even for equal masses seems to 
indicate that the fast sound and the concentration dependence of the gap are 
not related. 
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